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Invariance under local unitary operations is a fundamental property that must be obeyed by every 
proper measure of quantum entanglement. However, this is not the only aspect of entanglement 
theory where local unitaries play a relevant role. In the present work we show that the application of 
suitable local unitary operations defines a family of bipartite entanglement monotones, collectively 
referred to as "mirror entanglement". They are constructed by first considering the (squared) 
Hilbert- Schmidt distance of the state from the set of states obtained by applying to it a given 
local unitary. To the action of each different local unitary there corresponds a different distance. 
We then minimize these distances over the sets of local unitaries with different spectra, obtaining 
an entire family of different entanglement monotones. We show that these mirror entanglement 
monotones are organized in a hierarchical structure, and we establish the conditions that need to be 
imposed on the spectrum of a local unitary for the associated mirror entanglement to be faithful, 
i.e. to vanish on and only on separable pure states. We analyze in detail the properties of one 
particularly relevant member of the family, the "stellar mirror entanglement" associated to traceless 
local unitaries with nondegenerate spectrum and equispaced eigenvalues in the complex plane. This 
particular measure generalizes the original analysis of [Giampaolo and Illuminati, Phys. Rev. A 76, 
042301 (2007)], valid for qubits and qutrits. We prove that the stellar entanglement is a faithful 
bipartite entanglement monotone in any dimension, and that it is bounded from below by a function 
proportional to the linear entropy and from above by the linear entropy itself, coinciding with it in 
two- and three-dimensional spaces. 

PACS numbers: 03.67.Mn, 03.65.Ud 



I. INTRODUCTION 

Achieving a satisfactory understanding of the nature 
and structure of quantum correlations is ofparamount 
importance in quantum information theory (l| as well as 
in the study of complex quantum systems |2j; very re- 
cent speculations even hint at possible fundamental roles 
played by quantum entanglement in biological systems 
and processes with enhanced properties of quantum co- 
herence 0]. For bipartite quantum systems prepared in 
globally pure states, there is universal consensus on the 
fact that quantum correlations identify with entangle- 
ment, and they can be signaled by either one of several 
features that distinguish them from the classical ones. 
For instance, their non-local character, or their opera- 
tional significance in quantum informational primitives 
such as entanglement creation and distillation, or the per- 
formance they enable as resources for quantum communi- 
cation protocols such as teleportation [l[ . Entanglement 
in this case is simply measured by how much informa- 
tion one is missing on the state of the global system by 
accessing only one part of it, hence capturing the cor- 
relations between the two parties. Mathematically, such 
an information is quantified by the entropy of entangle- 
ment, i.e. the von Neumann entropy (or any monoton- 
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ically increasing function of it) of the reduced density 
matrix of either of the two subsystems [3], although it 
is worth remarking that other inequivalent measures of 
pure-state bipartite entanglement can be introduced be- 
yond the von Neumann entropy, such as the infinite set 
of Renyi entanglement entropies of the reduced density 
matrix (entanglement spectrum) [f| and the bipartite ge- 
ometric entanglement [fj , that are particularly relevant in 
the investigation of ground-state properties of condensed 
matter systems and in the theory of quantum information 
with continuous variables Q- 

For bipartite mixed states and for general multipartite 
states (pure or mixed), several inequivalent entanglement 
monotones (as well as a number of measures of more gen- 
eral types of quantum correlations beyond entanglement 
in mixed states Q) have been proposed, each apt to cap- 
ture a different signature of quantumness and/or possess- 
ing a different operational meaning, especially in relation 
to different informational tasks Whilst all entangle- 
ment monotones must vanish on separable states, there 
can be some that vanish also on some entangled states, if 
the latter fail to encode the particular resource character 
associated to a given entanglement measure. This is the 
case, for example, of bound entangled states, which have 
a nonzero entanglement cost but a vanishing distillablc 
entanglement, as no entangled singlets can be extracted 
from them by local operations and classical communica- 
tion (LOCC) in the asymptotic regime |Tc|. Fundamental 
requirements for any bona fide entanglement monotone 
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are then monotonicity under LOCC - operations that 
cannot increase bipartite entanglement on average - and 
invariance under local unitary (LU) operations pj| . The 
latter is in fact a basic requirement for any measure of 
correlation in general, since the choice of the local basis in 
which the density matrix of a system is expressed cannot 
of course affect the information encoded in shared cor- 
relations between two (or more) subsystems. This fact 
has motivated the search for the simplest LU-invariant 
"normal forms" of quantum states in discrete as well as 
continuous variable systems [l2j , so as to reduce as much 
as possible the number of state parameters needed for a 
complete evaluation of a particular entanglement mono- 
tone. 

Invariance of entanglement under LUs has however 
more far-reaching consequences, especially, and not en- 
tirely without surprise, in its quantification. Given a 
generic state (pure or mixed) of a bipartite systems, Li- 
Bin Fu investigated the consequences of applying a local 
cyclic operation on one of the subsystems [l3(. Fu de- 
noted by local cyclic operation any LU that leaves the 
corresponding reduced state invariant. Consider a bipar- 
tite quantum system (A\B) in a global pure state \iff) 
with reduced density matrices pa and pb respectively 
for subsystems A and B, and denote by Ua a LU acting 
on pa only. Then, the LU U a is cyclic if [Ua,Pa] = 0. 
Although being local, thus leaving the entanglement un- 
changed, such an operation changes the global state, 
yielding a nonlocal effect that can be detected only by 
measuring the two subsystems jointly. Employing the 
fidelity induced by the Hilbert-Schmidt metric, Fu iden- 
tified such a nonlocal effect by the distance between the 
initial and the final state. Recently, Fu's pioneering work 
has been greatly extended by Gharibian, Kampermann, 
and Bruss 

[3. They focused on the maximization of 
the Fu distance for bipartite states in Hilbcrt spaces 
of arbitrary dimension as a possible indicator of non- 
local properties, and derived and discussed closed formu- 
lae for the maximal Fu distance in three relevant cases: 
(pseudo)pure quantum states, Werner states, and two- 
qubit states. In between, two of us (S.M.G. and F.I.) 
investigated independently the consequences of LUs on 
global pure states of 2 x D and 3 x D bipartite sys- 
tems [ljf. They attacked the problem from the oppo- 
site side, investigating the minimization of the (squared) 
Fu distance, and proved, somewhat surprisingly, that in 
these two particular cases the minimum (squared) Fu dis- 
tance is a full bipartite entanglement monotone, coincid- 
ing with the linear entropy of entanglement (also known 
as tangle for qubit systems [H, Q3 ) ■ A similar analysis 
was performed to define a geometric LU-based entangle- 
ment measure for Gaussian states of continuous variable 
systems, where subsystem A comprises a single bosonic 
mode HI]. 

In the present paper we will present a complete gener- 
alization of the analysis carried out in Ref. |15( to all pure 
states of bipartite quantum systems with Hilbert space 
Hab of arbitrary finite dimension. We will show that one 



can construct an entire family of bipartite entanglement 
monotones that capture quantum correlations as quanti- 
fied by the action of minimally perturbing LUs on global 
pure states. Specifically, we prove that the (squared) 
Hilbert-Schmidt distance between a pure bipartite state 
IV- 1 ) ab G Wab = Ha ® Hb and the pure state Ua^ab 
obtained from it by applying a LU operation Ua on sub- 
system A only, once suitably optimized (minimized) over 
all LUs with fixed spectrum, defines a hierarchy of bipar- 
tite entanglement monotones. Moreover, the cyclic con- 
dition [Ua, Pa] = 0, rather than being imposed a priori is 
derived as a consequence of the minimization (optimiza- 
tion) procedure. We denote any such pure-state bipar- 
tite entanglement monotone as "mirror entanglement". 
This is pictorially reminiscent of someone mirroring her- 
self/himself in a mirror (which is curved, symbolizing the 
action of a LU): in absence of entanglement, the mirror 
image of the original pure state under a LU (the mirror) 
is a perfect reflection. Viceversa, the more entanglement 
is contained in the state of the system, the more distorted 
is the image that is reflected by the mirror. 

As just stated, no constraints on the employed LUs 
need to be imposed ab initio: all the mirror entanglement 
measures are proven to be LOCC-monotoncs for arbi- 
trary pure states of bipartite systems in any dimension. 
However, we will show that by imposing, as originally 
done in Ref. [l5[, specific requirements on Ua, namely 
a fully nondegenerate spectrum, one restricts the class 
further to faithful mirror entanglement monotones that 
vanish if and only if the state \i/j}ab is a product state. 
Upon further restricting the admissible LUs by requiring 
that they be traceless and with equispaced, nondegener- 
ate eigenvalues (thus with a pattern resembling a star in 
the complex plane), we will single out a special faithful 
mirror entanglement monotone, that we will name "stel- 
lar mirror entanglement" . We will prove that the stellar 
mirror entanglement enjoys the property of being a lower 
bound to the linear entropy in any dimension, reducing 
to the latter in the special cases Ha = C 2 and Ha = C 3 
originally considered in (l5| . Moreover, it is an upper 
bound to a function proportional to the linear entropy, 
with the proportionality constant being a simple function 
of the dimension of the reduced Hilbert space Ha- 

We remark that, by construction, the class of mirror 
entanglement measures is experimentally accessible by 
means of interferometric schemes fl9j involving at least 
two copies of a given bipartite state, one of which needs to 
be rotated by suitable LUs. In principle, each mirror en- 
tanglement monotone can be straightforwardly extended 
to mixed bipartite states via the conventional convex roof 
construction. Solving the convex roof optimization prob- 
lem is of course in general a formidable task. However, 
as the mirror and stellar entanglement measures are de- 
fined in terms of distances (in partial analogy with the 
construction of the geometric measures of entanglement 
H3])i m t ne conclusions we will briefly discuss how 
it might be possible to envisage alternative strategies to 
compute their mixed-state extension without resorting 
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directly to the convex roof construction. 

The paper is organized as follows. In Scc.|n]wc define 
the class of mirror entanglement measures for all pure 
states of bipartite quantum systems, prove their mono- 
tonicity under LOCC, characterize their hierarchic struc- 
ture, and determine the conditions for faithfulness. In 
Sec. Mil we focus on the stellar mirror entanglement and 
investigate its relationship with the linear entropy of en- 
tanglement, providing the exact lower and upper bounds 
that relate the two quantities (The detailed proofs of 
the bounds are reported in two Appendices). Finally, in 
Sec. II VI we briefly discuss some of the implications of our 
results and some possible future research directions con- 
cerning the extension to mixed states and the problem 
of identifying total and partial pure-state factorization 
(separability) in quantum many-body systems. 



II. MIRROR ENTANGLEMENT: DEFINITION 
AND MONOTONICITY 



Let us consider a bipartite system in a pure quantum 
state = \tP)ab belonging to a Hilbert space Hab = 
Ha® Kb = < & dA ® C dfi . We will assume without loss 
of generality that d = d,A < ds- Let us consider general 
LUs acting on Ha of the form 

W K a = Wk = Y1 x MM^\, (1) 



where 



A. = {\ j } = {e ie >} (j = !,...,<*), (2) 



denotes the spectrum of the eigenvalues of W A . The max- 
imal fidelity (squared overlap) between state l^) and the 
LU-transformed state WaI^) is 



maxKV>|W A M| s 



Wa 



(3) 



Definition ([4]) is recast in terms of the reduced state (JH) 
by rewriting Eq. ^ as 

F^ = max|Tr[W A |V')(V'|]| 2 
= max|Tr A [VF A Tr B 

W A 

= max|Tr[W A e]| 2 . 
W A 



By the monotonicity of the square function one can write 

J F§ = max^ A |Tr[W A £]|. Let \i) be the eigenbasis of 

g and pi its eigenvalues, so that one has the spectral 
decomposition g = The set of allowed LUs, 

that is, the set of unitary matrices acting on Ha with 
spectrum A [Eq. (j2l], can be written in terms of Va — 
£\ Xi\i){i\ as 



(6) 



where U rotates the eigenbasis of g into the eigenbasis 
of Wa, \4>i) = U\i). In principle, U can be any SU{d) 

unitary matrix. We can write 



max |Tr \UV\U* q] I 

U^SU(d) 1 1 



max 

uesu(d) 



max 

uesu(d) 



E^E 



Pj\Uij\ 



(7) 



and it takes non- negative real values in the interval [0, 1]. one nas | ^ ^ p. 



where Uij = (i\U\j). 

We will now show (Theorems [5] and 0] below) that the 
AME measures are indeed legitimate pure-state entan- 
glement monotones. 

Theorem 2 The ME vanishes on pure separable (i.e., 
product) bipartite states \ip®) = \iPa) ® \^Pb)- 

Proof. For a product state, g is a rank-one matrix, with 
eigenvalues pj = Sjk for some index fc. Choosing [7 = 1 

lA fc | = 



Si Ai V , (\,t,(),i | 



Definition 1 (mirror Entanglement) The bipartite 
A-mirror entanglement (AME) between subsystems A 
and B in the state \tp) is defined as the square of the 
minimum Euclidean distance between \ip) and the set of 
transformed states Wa\iP) obtained by the action of LUs 
of the form Eq. (QP with spectrum A on subsystem A: 



e l6k \ = 1 = JF A . Therefore, from Eq. Q £ A (>®) = 



£ A W0 =min(l- MWaW 

Wa 



1-^ 



(4) 



Consider the reduced density matrix of subsystem A: 

Q=QA = ^ B Mm ■ (5) 



for any A. □ 

Before tackling the monotonicity of the ME under 
LOCC, we first prove an auxiliary lemma that simpli- 
fies the optimization problem involved in the definition 
of the ME. 

Lemma 3 The maximizing unitary U in Eq. ^ is a 
permutation matrix. 

Proof. Eq. ([7]) can be written as 

! F A = max \Tr\M$B(U)}\ , 
U£SU(d) 1 v 1 
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where [B(U)]i 



I 2 and (M#) 



PiXj. Noticing 



that B(U) is a unistochastic matrix, we can write 



>Ff= max \Tr[M$B}\ < max \Tr[M^B] 

v B unistoch. 1 BeB d 



where we have enlarged the optimization domain to the 
whole set of all d x d doubly stochastic matrices, 
i.e. the d-dimensional Birkhoff polytopc. By the Birkhoff- 
von Neumann theorem, is the convex hull of the set 
Sd of d x d permutation matrices (that is, the permuta- 
tion matrices in Sd are the extreme points of B f /). We 
can thus write B = J2k QkSk, where Sk € Sd, and q = 
{qk} is a d!-dimensional probability vector. The maxi- 



mal fidelity becomes 

max f Efe% T±[M$S k ] 
angle inequality. Let S, 
trix that maximizes 



< 



Ek1kTr[M$S k ] 

where we have used the tri- 
be the permutation ma- 
Then , [fJ < 



T±[M$S a 



We are 



(max^ fc %) Tr[M^ ma 

left to show that Smax = B(U) for some U. This is 
achieved by noticing that all permutation matrices, in- 
cluding S'max, are orthogonal and hence unitary, and that 
B(S max ) = 5 max . This concludes the proof. □ 

As a corollary of Lemma [3J we find that the optimal 
LU operation W\ = S max V\S^ nax that maximizes 
[Eq. ©] commutes with the reduced state g. To see this, 
let us write S'max = J2i l CT i)(*l> where a is the permutation 
described by the matrix S max . Then 

Wa = ^><><i|^ Va (z>'>fo 

= ^A i ki)(a i | (8) 

i 



which is diagonal in the same basis as g, and therefore 
[g, W\) = 0. The last result shows that the eigenvectors 
of the optimal LU Wa that solves the minimization in 
the definition of the AME [Eq. ((H)] are just obtained by 
a reordering of the eigenvectors of the reduced state g. 

Collecting the previous findings, \ F§ can be written as 



F£ = max |Tr[MAs] 



ses d 



(9) 



We are now ready to prove the following important 
result. 



Proof. We will prove 1 that F^ is monotonically in- 
creasing under LOCC, which implies the statement. An 
arbitrary pure state \ip) is ensemble-transformed under 
LOCC according to: \ip) — > {pi, where each LOCC- 

transformed state of the ensemble reads: 



/p^ipi) = (A, ® l B )\ip) 



(10) 



The positive weights {pi} satisfy the normalization con- 
dition pi = 1, while the Kraus operators associated to 
the local dynamics satisfy the POVM (positive operator 
valued measure) completeness relation: ^ t A\Ai = 1. 
Let the reduced state be p = Ttb[\iP) (tp\] and likewise the 
reduced LOCC-transformed states be p { = Tr B [|V>i)(V'i|]- 
For each of the latter, the local dynamics yields: 



Pi Pi 



A iP A\ 



(11) 



Then, in order to prove that £(i/j) > YliPi^i^iJi it is 
sufficient to prove that . pi y/F(ipi) > -J F(ip). Let 



M\fp = V AipAjVi = Jpisfp^i 



(12) 



be the polar decomposition of Ai^fp, where Vj is a suit- 
able unitary matrix. Exploiting the properties of F, we 
can write: 



Pi vf(^) = J2 pi ^ x i^iPiWi] i 

i i 

= J2p i nmx\Tr[p i ViW i V? 

i 

= ^p i ma X |Tr[(^ t Vp i )(Vp i KO^] 

i 
i 

>max^|Tr[Vp4^VpW] 

i 



> max 
w 



>xnax\Ti:[pW}\ 



(13) 



This concludes the proof that local operations on party 
A do not increase 1 — F. The proof for operations on 



Theorem 4 The ME is monotonically non-increasing 
under LOCC operations, i.e., is a full pure- state bipartite 
entanglement monotone. 



1 Here and in the following steps of the proof we omit the label A 
for ease of notation. 
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party B follows trivially. We have: 

y^Piy/Ftyi) = VV max|(^|T^(g)lL|V>i) 

— — vVi 



> max 



(14) 



Therefore, since £ (ip) = 1 — and all local operations 
do not increase 1 — F, we have proved that the ME is 
non-increasing under LOCC. □ 

We have thus introduced a family of pure-state bipar- 
tite entanglement monotones, that satisfy the three fun- 
damental axiomatic properties of vanishing on separable 
states, being invariant under local unitaries, and being 
monotonic under LOCC fl ©, GH, H HI . They are as- 
sociated to the minimum distance between a quantum 
state \rp) and its image after suitable unitary operations 
with fixed spectrum performed on one subsystem only. 
Surprisingly enough, one sees that a LU operation on one 
part of a bipartite system, while leaving the entanglement 
invariant, leads nevertheless to a state alteration whose 
proper quantification is itself an entanglement measure. 
The properties of the spectrum A define the shape, or re- 
flectivity, of a fictitious mirror which produces the image 
of the quantum state after the action of a LU. In the ab- 
sence of entanglement, there exists always one such LU 
that leaves the state invariant, yielding a perfect reflec- 
tion from the mirror. If the state is entangled, the 
action of the minimal or least perturbing LU with spec- 
trum A necessarily results in a transformed state which 
has a nonmaximal fidelity with in turn, this dis- 
tortion quantifies the amount of bipartite entanglement 
present in \tp). 

The class of AME exhibits a hierarchical structure 
depending on the characteristics of the spectrum A = 
{e z9j }. One extreme case is represented by 6j = Vj. 
In this case the identity is clearly the extremal LU op- 
eration that defines the 1ME according to Eq. ((4]), and 
the ensuing entanglement measure £t(ip) is trivially zero 
for all quantum states |?/>). Progressing in the hierarchy, 
if A contains an r-dcgcncratc eigenvalue, then the corre- 
sponding ME vanishes on entangled states \tp) of Schmidt 
rank < r. On the opposite extreme, if the eigenvalues in 
A are all nondegenerate, one obtains the most sensitive 
measure of mirror entanglement AME. This classification 
is summarized by the following theorem: 

Theorem 5 Let Wa be a LU and let the associated spec- 
trum A have degeneracy r ( any eigenvalue repeated at 
most r times). Then £a(V>) = if and only if the Schmidt 
rank of SR(?/>), is no larger than r. 

Proof. The Schmidt rank of \ip) amounts to the rank 
of the reduced density matrix g, or the number of non- 



vanishing elements in the probability vector p. We first 
prove sufficiency and then necessity: 

I) [SR#) < r => £ A (V0 = 0]. 

If the Schmidt rank of ip is s < r then 
|Tr[S , VA<5'^]| = 1 can be attained by any per- 
mutation matrix 5 max such that S^ aK maps the 
s-dimensional domain of g into the r-fold degen- 
erate subspace of V\, so that |Tr[S'VAS'^gi]| = 
\J2i=i ^iPi\ = 1j where Ai = ... = A r . In this 
way one has £a(V0 = 0. 

II) [£ A (i>) = SR(V0 < r}. 

To prove the inverse implication, let a be the max- 
imizing permutation in Eq. Aj = A -i, and let 
E be the set of indices for which pi ^ 0. Then 
s 



|E|. Expanding F§ we get 



< 



Erf 
Erf 
Erf 



Ai Pi 



E X i X oPiPj 

{*W}e£ 



53 (AJAj- + \j\i)piPj 
{i<i}es 

2 E P* p i 
{i<i}es 

2 

= 1 . 



The inequality follows from (a*b + b*a) < 2|a||6|, 
and it is saturated if and only if Re[a*fe] = 1, thus 
for each pair i, j 6 S we have 



A.; = Xj \fi, j G E 



(16) 



By assumption, there are at most r equal values of 
A, thus |E| < r. This completes the proof. □ 

By proper parameterization one can obtain a picto- 
rial representation of all the AME monotones. To do so, 
consider again the LU spectrum A = {e l6i }. The AME 
is clearly invariant under permutations of the phases Qi , 
therefore one can restrict to < 6\ < 62 < . . . < £>d < 27r 
without loss of generality. Furthermore, every monotone 
£a is invariant under global phase shifts of the eigenval- 
ues, thus 6\ = can be assumed throughout. Finally, 
one can equivalently re-parameterize the phases by the 
set of gaps ip k = (Ok+i - k )/2Tr, with ip d = I - 9 d /2ir. 
Since the phases are ordered, one has tp k > and clearly 
J2k=i 'Pk — 1 [sec Fig. [I]. Thus, each AME is in corre- 
spondence with a d-dimcnsional probability distribution 
{(fik}- This means that all the AME monotones can be 
represented in a simplex, as shown in Fig. [2] 

As a Corollary of Theorem[5j one can easily identify all 
the AME monotones which arc faithful, i.e., vanish on and 
only on product states. These correspond to fully non- 
degenerate spectra A, and fill the entire shaded region in 
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FIG. 1: Parameterization of the spectrum A of the LU Wa, 
given by the normalized phases 9i [Eq. ([2j|] or the gaps <fi 
between phases (in units of 27rrads). A spectrum is degen- 
erate whenever a ifk equals zero, which corresponds to the 
boundary of the simplex [see Fig. [2] . 



Fig. [2] A particular case of a faithful measure is obtained 
when the eigenvalues constituting the LU spectrum A are 
cquispaccd in the complex plane and add up to zero, as 
it is the case, e.g., when they are taken to be the dth 
roots of (— 1) . Since their representation in the Ar- 
gand diagram resembles a star (or a regular polygon), 
we name the associated entanglement monotone "stellar 
mirror entanglement" (-ArME), or "stellar entanglement" 
for brevity. The properties of this particularly important 
monotone are discussed in the forthcoming Section Hill 

For every mixed state qab of a bipartite quantum sys- 
tem, and for any A, the AME can be defined via the 
convex roof construction 



Sa(qab) = inf PiSA(ipi) 
{Pi,i>ii 



(17) 



where the infimum runs over all pure-state decomposi- 
tions of qab = '^2iPi\ipi)('>Pi\- The convex roof construc- 
tion ensures, by definition, that the ME is still a full 
LOCC monotone, and hence a valid measure of entan- 
glement, even in the general case of mixed states. It 
is obvious that, like in the case of most other measures 
of entanglement, the actual computation of the convex 
roof on generic mixed quantum states is a formidable 
task. In the conclusions we will briefly discuss how the 
ME measures, being defined in terms of geometric dis- 
tances in Hilbert space, might be extended to the general 
mixed-state case via alternative procedures that bypass 
and avoid the explicit evaluation of the convex roof con- 
struction. 



FIG. 2: (Color online) Representation of the set of A-mirror 
entanglement monotones for d — 3 in the three-dimensional 
simplex of the gaps {ipk} associated to the spectrum of A. 
The stellar monotone £^ is placed at the center of the sim- 
plex (ipi — 1/3) and denoted by a full (red) circle. The 
trivial monotones, i.e. those that are identically vanishing, 
are associated to spectra for which two out of three ipi's are 
zero. These monotones fall at the extremal points (corners) of 
the simplex and are represented by full (white) circles (white 
dots). Partially degenerate monotones £a are those for which 
one or more ipt is vanishing. They fall at the boundaries of the 
simplex. For instance, the particular case reported by a full 
(green) circle on one side of the simplex represents a partially 
degenerate monotone that is nonvanishing only on entangled 
states with Schmidt rank SR=3. 



Section, characterized by a LU with spectrum A [Eq. ©] 
obeying precise constraints, defined in the following. Ob- 
serving that the LUs entering in the optimization for the 
AME can always be represented as W\ = UV\W , as in 
Eq. (jnj), we then define the stellar spectrum A = it as 
the one such that 



exp 



a 



(18) 



where the matrix S z represents the z component of the 
spin- J operator with J = (d — l)/2, S z = diag{ J, J — 
1,...,— (J — 1),— J}. Explicitly, the stellar spectrum is 
given by the diagonal of , 

* = {e"'}, 0j = d ~ 2] d +l * (j = l,.-.,d). (19) 

Definition 6 (Stellar mirror Entanglement) The 

stellar mirror entanglement (~kME) £^(tJj) is defined as 
the AME [see Definition^ Eq. ^] with A = ~k . 



III. STELLAR MIRROR ENTANGLEMENT 

In this Section we study in more detail a special in- 
stance of the AME monotones introduced in the previous 



Reminding that the local dimension d is defined as 
d = minld^jds}, the eigenvalues in the stellar spec- 
trum correspond to the dih complex roots of (— 1) x , 
do not give rise to any degeneracy, and are equispaccd 
in the complex plane (resembling a star once connected 
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by rays to the origin). It is straightforward to show that 
Ylj=i e%9i = 0; as a consequence, all the LUs W^. with 
spectrum ir are traceless: TrfW^] = 0. According to 
Theorem [5] proved above, the stellar mirror entanglement 
(or, for brevity, stellar entanglement) £^(ip) vanishes if 
and only if \ijj) is a separable (product) state; thus, the 
■ArME is a faithful entanglement monotone. This result 
includes and extends to any dimension d the correspond- 
ing finding of Rcf. [15j . valid for d — 2, 3. It also gener- 
alizes to any dimension the LU-based separability crite- 
rion that was established in that same paper for d = 2, 3. 
Indeed, consider the optimal LU , i.e. the Single- 
QUDit Unitary Operation (SQUDUO) that realizes the 
stellar entanglement by minimizing the (squared) Eu- 
clidean distance over the entire set of LUs with stellar 
spectrum: £^{ip) = minimi (l — |(V'|W-A-IV')| 2 )- ft i s 
straightforward to show that the faithfulness of the stel- 
lar entanglement £^.(ip) implies the following LU-based 
separability criterion: A pure state \ip) of a bipartite 
system is separable (product) if and only if the optimal 
SQUDUO W7* leaves it invariant: W^\ip) = \ip). The 
faithfulness of what we here call TirME, when restricted to 
the qubit case d = 2 , has played a key role in the de- 
velopment of a general theory for the exact and rigorous 
detection and characterization of fully factorized ground 
states in several classes of non-exactly solvable spin-i 
models on translationally-invariant lattices [23[ as well 
as in more general geometries [24| with arbitrary spatial 
dimensions, both at finite size and in the thermodynamic 
limit. Based on the general results of the present work, 
in the concluding Section IIVI we will discuss some possi- 
ble guidelines for the extensions of such methods to the 
problem of the occurrence of total and partial factoriza- 
tions (such as dimerization, trimerization, and polymer- 
ization) in models with local spin variables of arbitrary 
dimension. 

At this stage, we notice that by using the results of 
Lemma[31 according to which the optimal change-of-basis 
matrix U is a permutation matrix, we can write the fol- 
lowing compact expression for the -frME: 



£^{ip) — min 1 



d 



(20) 



where the optimization is over all permutations a encod- 
ing a reordering of the eigenvalues pk of the reduced state 
g. Equipped with this expression, we can proceed to in- 
vestigate the relation of the ■rfrME with other measures of 
entanglement that are expressed as sums of products of 
eigenvalues of the reduced density matrix. The foremost 
measure of this kind is the linear entropy of entangle- 
ment E L {il)), defined as the linear entropy of the reduced 
density matrix g (22|: 



E L ty) =S L {g) 



1 ^ 
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FIG. 3: (Color online) Behavior of the stellar entanglement 
£^ plotted against the linear entropy of entanglement El for 
20000 random pure bipartite states ip £ <C dA g) C dB with local 
dimension d — min{dA, ds} = 4. The upper boundary (black 
online) is given by the bisectrix £ ^ = El ■ These states, that 
accommodate for a stellar entanglement coinciding with the 
linear entropy are characterized by a three-fold degenerate 
eigenvalue in the spectrum of the reduced state g: p\ = P2 — 
P3 = I, pa = 1 — p (with < p < 1). The lower boundary 
is branched into two different segments. For < El < 2/3, 
the minimum £^ satisfies £^ = \El\ this boundary (red 
online) accommodates states \%p) with rank-2 marginal states 
q: pi = p, p2 — 1 — p, P3 = Pa = 0. The second branch 
(blue online) of the lower boundary, defined for 2/3 < El < 
1, accommodates states satisfying £^ = \El — \- These 
states have reduced density matrix q with doubly degenerate 
spectrum of the form: p\ = p2 = i^ 2 , P3 — pa = All 
the quantities plotted are dimensionless. 



(21) 



It has been observed in Ref. [15| that the -frME and 
the linear entropy of entanglement coincide exactly on 
all pure states \ip) of bipartite systems with reduced 
states g of local dimension d — 2 or d = 3 (qubit 
or qutrit). This can be easily verified by comparing 
Eqs. (|2U)) and (f2"Tj) and recalling the normalization con- 
dition Tig = X)iP« = 1. In general, however, the two 
measures can differ, resulting in an inequivalcnt ordering 
imposed on the set of pure entangled quantum states \tp) 
with respect to a bipartition involving at least a qudit 
with d > 4. It is interesting to illustrate explicitly the 
discrepancy between the linear entropy and the stellar en- 
tanglement in the case- d = 4. In Fig. [3] we report £Ll ver- 
sus El for a sample of 20000 randomly generated states 
ip 6 C dA <g) C dB with d A — 4 and arbitrary d B > 4 (upon 
applying the Schmidt decomposition, the effective dimen- 
sion of each subsystem is reduced to d = min{ci J 4, ds} 
which amounts to 4 in this example). We find that, al- 
though physical qurfit states span a two-dimensional sur- 
face in the (El,£^) plane, nonetheless for a given values 
of El there exist sharp upper and lower bounds on £^. 
In fact, we will see below that these bounds are general 
and admit an exact analytical expression in any dimen- 
sion. The classes of entangled states that saturate them 
in the case d = 4 are specified in the caption of Fig. [3] 
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The pattern observed for the local dimension d < 4 
extends indeed to arbitrary values of the local dimension 
d. In particular, the hierarchical relationship £^ < El 
always holds (with equality for local dimension d = 2,3), 
and a rigorous proof is provided below. Moreover, the 
stellar entanglement icMF, presents a structured, multi- 
branched lower bound as a function of the linear entropy 
of entanglement, with the number of branches growing 
with local dimension d > 4. Without aiming at a char- 
acterization of the complete lower boundary for states in 
Hilbert spaces of arbitrary dimension, we can neverthe- 
less show that a simple rescaling of the linear entropy 
El allows to derive an exact analytical lower bound on 
£^ that holds in general for arbitrary dimension. In the 
particular case of local dimension d = 4 it corresponds to 
the bottommost branch of the lower boundary in Fig. [3] 
Indeed, one can prove that the following holds: 

Theorem 7 The ~kME is a lower bound to the linear en- 
tropy of entanglement and is an upper bound to a rescal- 
ing of it on all pure bipartite states in Hilbert spaces of 
arbitrary finite dimension: 



2{d- l)sin 2 (7r/d) 



E L ty) <£+{i>) <E L m . (22) 



The rightmost inequality is always tight and the states 

that saturate it are those for which | \Pi\, with 
being the eigenvalues of the stellar spectrum, Eq. (JTHJ) , is 
invariant under permutations. Incidentally, we note that 
for local dimension d = 2,3, all pure states have such 
permutational invariance. The proof of the rightmost 
inequality is provided in the Appendix [X] 

The leftmost inequality is generally tight for El < 
d/[2(d — 1)] and is saturated by states \ip) with rank- 



2 marginal states g: p a 



P, Po 



1 ~P, Po 







[i = 3, ...,d), for any local dimension d. For d = 2,3, 
the lower and upper bounds in Eq. (|22[) coincide. The 
leftmost inequality is proven in Appendix [Bj 

The attainable region for £^ at a fixed El increases 
with increasing local dimension d, and the lower bound 
vanishes asymptotically as d — ¥ oo, showing that there 
can exist pure bipartite entangled states of two qucfits 
with local dimension d 3> 1, whose linear entropy of en- 
tanglement lies in the range < El < \ 1 and yet pos- 
sess an infinitesimal degree of stellar entanglement -frME. 
The situation may change again and the trend may be 
reverted in the infinite-dimensional case. Indeed, in a 
previous study, some of us have shown that restricting to 
Gaussian pure bipartite states of two continuous variable 
modes, and to Gaussian LU operations, one can define a 
specific Gaussian counterpart to the stellar entanglement 
TirME that amounts to a simple monotonically increas- 
ing function of the linear entropy of entanglement fl8j . 
and thus back in analogy to the case of low-dimensional 
qudits (d = 2,3) [l^]. In conclusion, the analysis re- 
ported in this Section shows that the stellar entanglement 
tVME, a prominent representative of the AME mirror en- 



tanglement monotones, provides in general an indepen- 
dent characterization of bipartite entanglement in pure 
states of arbitrary finite dimension and is endowed with 
an intrinsic geometric origin depending entirely upon the 
global non-local effects induced only by the action of suit- 
ably identified classes of LU operations. 



IV. DISCUSSION AND CONCLUSION 

In the present work we have introduce a geometric 
framework to derive bipartite entanglement monotones, 
including faithful ones, that apply to all pure states in 
Hilbert spaces of arbitrary dimension. The measures of 
mirror entanglement and the faithful mirror stellar en- 
tanglement that we have introduced are defined in terms 
of the minimal (squared) distance from a pure state to 
the pure state obtained from it by the action of suit- 
ably optimized LUs acting only on one part of the bipar- 
tite quantum system under consideration. We identified 
a hierarchy of these LU-based entanglement monotones, 
studied their properties, and determined conditions for 
their faithfulness. Among the faithful mirror entangle- 
ment measures, we focused on a special instance, the 
stellar mirror entanglement, characterized by additional 
symmetries in the spectrum of the associated LUs. We 
proved that the stellar mirror entanglement obeys upper 
and lower bounds as a function of the linear entropy of 
entanglement in any dimension, reducing to the latter for 
local dimension d = 2,3. Our results generalize an earlier 
study limited to pure states with reduced density matri- 
ces of lower local dimension d = 2, 3 (l5j. Our work goes 
along a complementary direction compared to other stud- 
ies that have investigated the nature of non-local effects 
when the distance from a state and its image under the 
action of LUs is maximized rather than minimized [TJ] . 
It is remarkable and somewhat surprising that looking at 
such a simple structure of LUs from opposite ends can 
provide so much insight both on the structure of pure 
state bipartite entanglement and, at the same time, on 
the patterns of non-local effects generated from the appli- 
cation of LU operations. This interplay, yielding a host of 
complementary results, might hide yet undiscovered fea- 
tures common to the two approaches. Both rely on the 
natural physical intuition of the operational approach to 
the study of physical systems: looking at the response to 
a given action as basic diagnostic tool of physical prop- 
erties. It might then be worth to try to compare the two 
different situations in terms of a classification/ordering 
of least-disturbing and/or maximally disturbing LUs. 

An interesting and important subject for future re- 
search is the generalization of the entanglement structure 
generated by LUs to include mixed states. It has been 
recently shown [25[ that the problem of computing the 
convex roof of a prototype distance measure of entan- 
glement such as the global geometric entanglement [2(| 
can in fact be recast in terms of determining the maxi- 
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mal fidelity to a separable state. It is tempting to specu- 
late that this important result might perhaps be adapted 
and generalized to other classes of distance measures of 
entanglement, such as the mirror entanglement and the 
stellar entanglement introduced in the present work, or 
the mixed bipartite-multipartite geometric measures of 
entanglement defined as the hierarchy of distances from 
the sets of inseparable states . 

Extensions of the present investigation to the qualifi- 
cation and quantification of multipartite entanglement 
and the characterization of monogamy constraints on 
its distribution appear to be challenging, even when re- 
stricted to pure states. Here the problem appears of 
course to be that of understanding the nature of gen- 
eralized "local operations" in the multipartite setting, 
and their ordering according to the associated local di- 
mension. Viceversa, a more readily exploitable appli- 
cation of our results is concerned with the factorization 
of quantum ground states in cooperative qudit systems, 
a phenomenon that is currently receiving significant at- 
tention from both quantum information and condensed 
matter communities 0, [23|, [24|, [26|. Indeed, the variety 
of quantum states belonging to one and the same LU- 
cquivalence class may anyway have rather distinct fea- 
tures that become relevant especially in the context of 
many-body physics. For instance, if we consider a spin 
chain in a perfectly ferromagnetic state | tt ■ • • t) > m P~ 
ping every other spin amounts to a LU operation that 
has no effect on the entanglement, yet results in a totally 
different ordered phase with a vanishing magnetization 2 . 
The specific form of the ground state of a many-body 
Hamiltonian, then, is important as well as its entangle- 
ment content and distribution in the form of bipartite 
or multipartite quantum correlations Q. In this con- 
text, the formalism of LU-based geometric entanglement 
for states with reduced density matrices of lower local 
dimension d — 2,3 has already been applied to define 
energy witnesses as efficient diagnostic tools of factoriza- 
tion, relating it to the vanishing of faithful entanglement 
measures such as the stellar mirror entanglement. How- 
ever, this LU-based theory of ground-state factorization 
had to be restricted so far to spin- 1/2 systems [13] and 
to be limited to investigate only total factorization into 
products of single-spin states \2j^ [24j], as the extension 
to higher spin systems and partial factorizations (such 
as dimerization) required considering higher local dimen- 
sions d. The results of the present work allow in principle 
to extend the LU-based methods, originally developed 
for spin- 1/2 models and total factorization, to the study 
of total as well as partial factorization points in models 



2 In this context, we should also notice the fact that the states 
in the above example are fully factorized (fully separable), while 
the condensed matter terminology identifies them as strongly 
correlated, referring to the behavior of the two-point correlation 
functions; here we will adhere to the conventions of entanglement 
theory. 



of higher-dimensional spins and in spin- 1/2 models with 
frustration. In the latter case of frustrated spin-1/2 sys- 
tems, the LU-based entanglement formalism has been re- 
cently exploited to relate the existence of fully separable 
ground states (totally factorized in the tensor product 
of single-spin states) to the absence of frustration [28j |. 
Equipped with the general proof of equivalence between 
pure-state factorization and invariancc under LU spin op- 
erations with stellar (or, generally non-degenerate) spec- 
trum, one may now investigate the possible existence of 
dimerized quantum phases (i.e. ground states composed 
of factorized singlets: tensor products of d = 4 units, 
each of which is internally entangled) by looking for can- 
didate ground states invariant under stellar LUs in d = 4. 
For spin-1/2 chains in the thermodynamics limit or in 
2D lattices one may expect that a hierarchy of compat- 
ible types of ground state correlations may take place, 
ranging from full factorization into products of single- 
spin states, to dimerization, up to genuine multipartite- 
entangled phases, as the degree of frustration increases 
as a function of a given tunable external magnetic field. 
The exploration of these intriguing scenarios, as enabled 
by the general analysis developed in the present work, 
will be the subject of further future studies. 
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Appendix A: Proof of £^(4>) < El(iP) 

In order to prove the theorem, we will use Dirac bracket 
notation throughout. It is convenient to establish some 
definitions and facts. Let us denote the maximally mixed 
vector by |1) = (1/d, l/d,... , l/d)l With this, = 
1 jd and we denote the projector onto the 1 1) subspace by 
P 11 = ep)(l|, the complement being = 1 - pH. The 
eigenvalues A^ can be arranged in the vector |A), with 
complex conjugate |A*). We have (A|A) = (A*|A*) = d 
and for the stellar monotone, also (A|l) = 0. 

Define M a = ±<r T (|A)(A| + |A*><A*|)cr. Denoting the 
identity permutation by e, we can write M = M e = 
Re |A)(A|, and M a = a T Ma. It is clear that M > and 
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thus M a > for all a. Let 

g a {p) = 1 - I E X ^ I 3 = 1 - (P\ M -\P)- ( A1 ) 

i 

For any < s < 1, let £/ e ,s = {p|g e (p) > s}. Proceed now 
to define C s = {p \ £^{p) > s}. Then C s = r\ ae s d Q s ,a- 

Before proceeding any further, we will find the follow- 
ing relation useful, 



o"es d 



dl- 



d-1 



-P A 



(A2) 



which can be easily derived by use of Schur's lemma. 

The proof consists on proving a lower bound of El over 
the set C s . Equivalently, we may search for the minimum 
of the objective function f(p) = 1 — (p\p), 



minimize f(p) 
subject to p £ C s . 



(A3) 



Let us denote the solution to this problem by p* and the 
achieved value f(p*) by /* = f(p*). 

Lemma: The solution to the optimization in (|A3|) sat- 
isfies 



d-i 



(A4) 



Proof. Consider the vector |g) = (1 — y/1 — s)|l) + 
y/l — s|ei), where |ei) = (1,0,0,...)*. This vector has 
9<j{q) = s, for all a, thus q <E C s , and f(q) = s^p. 
Therefore f* < f(q) = s^=± □ 

Lemma: The solution to the optimization in (|A3|) sat- 
isfies 



d-i 



(A5) 



Proof. The program (|A3[) is spelled out as 

minimize / (p) 

[ s-g a (p)<0 
-Pi < 

$>-i = o' 



subject to 



(A6) 



and the Lagrangian is 

L(p,fx,a,v) = 
f(p) + Y^(x a (s- 9<T (p))-(a\p}-is(l-d(l\p)). (A7) 



Consider the dual program [29|,|30j to that of ([AB]). The 
dual Lagragian is 



L_d(/z, a, i/) = inf L(p, fi, a, v) 
p 



(A8) 



In order to have the dual Lagrangian Loiter, a, v) > — oo 
it is convenient to write the objective function as 



./•( P )H P |(d 2 |i)(i|-i)b). 

L(p, n, a, v) can then be written as 



(A9) 



L(p,n,a,v)=(p\ ^ 2 |1)(1| -l + ^ ff M ff j |p) 

(A10) 

-«a|-i/d(l|)|p>-(l 



Hence, L is bounded from below whenever H = V 2 L = 

fJ, a M a — 1 > 0. In such case, the infimum in 
can be readily computed by solving V p i = 0, 

1 



\p) = -H- 1 [\a)-vd\l)], 



(All) 



yielding 



M Ml a, v) = - - (H - vd(l\) H- 1 (\a) - ud\l)) 



(1 - s)^/i 



(A12) 



At this point, any value of Ld with {fi a } satisfying H > 
yields a lower bound to f*(s). Moreover, observe that 
J2 a ^=Tr[H]/d, 



L D (v, a,v)=-- ((a| - ud{t\)H- x (\a) - vd\t)) 



(1 - s)Tr[H]/d-u. 



(A13) 



We can choose to evaluate Lr>{yL^a,v) at \i a = [i so 
that H becomes 



(m^-i 



Defining x = fijzj — 1, we have 



(A14) 



H = xP 1 - + (d- 1)P 11 (A15a) 
H- 1 = x^P 1 - + (d- I)" 1 / 311 , (A15b) 

while the H > condition reduces to x > or /i > 
Moreover, = Tr[#]/d = (1 + x)<h£. Thus, we 

obtain 

L D {x,a,v) = - ^lll a± )ll 2 - C 1 ~ «)— j^a; 
1 



(1 -.)_-„. 



(A16) 
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Maximizing over x we obtain a;* = jjjjj^H *J ^rj an d 
L D {x*,a,v) = -111^)11^/(1-5)^ 



4(d-l) 



|a")-W|l)|| : 



-(1 -,)_-„. 



(A17) 



The term proportional to 1 1 1 ck -1- ^ 1 1 can only be negative, 
thus we set it to zero by choosing \a) — |a") = ay/d\t) 
(a > 0), leaving 

T r * \ (a-vy/d) 2 d-1 
L D (x*,a,v) = - V 4(d _ 1) ; - f - (1 - s )^-- ( Al8 ) 

Maximizing over v we get i/* = 2-^ — 4^p and 

L D (a;*, a ,i/*) = a *zi- " (A19) 
a y/d 

Since a > the maximum is achieved at a* = a = 
yielding 



Ld(x* ,a*,v*) = s 



d-1 



which shows that 



d-i 



Combining the last two lemmas, we have 

d-1 



a 

From this we have the following relation, 

S^ip) > s=>/(p) > s^—^, 
or equivalcntly 



(A20) 

(A21) 
□ 

(A22) 

(A23) 
(A24) 



E+(p) > s=>E L (p) > s. 

In particular, for points with £^{p) — s we have 

S+(p) = s < E L (p) . (A25) 

This completes the proof of the rightmost inequality in 
Theorem El □ 

Appendix B: Proof of 2^ sin 2 (Tr/d)E L (ip) < 

For this proof we will use the same notation as in Ap- 
pendix fAJ Additionally, let A = |A)(A|. The approach of 



the proof is the following. Since all entanglement mono- 
tones can be expressed as function of the spectrum of 
the reduced density matrix, we will talk about proba- 
bility vectors instead of quantum states. We will show 
that any probability vector \p) can be obtained by a 
series of transformations starting from the pure vector 
\q) = (1, 0, . . .)t such that the increments in both en- 
tanglement monotones under the action of these trans- 
formations always verify 

A£^>2^-j^sm 2 (iT/d)AE L . (Bl) 

This, combined with the fact that both monotones van- 
ish on (1, 0, 0, . . .) completes the proof. We now proceed 
to prove Eq. (|B1|) . An essential ingredient in the proof 
consists in obtaining a sequence of transformations that 
will bring q to p while still having a manageable form of 
Af^. To this end, we realize that the most disturbing 
element in the expression of are the optimizations 

over the set of permutation matrices for the initial and 
final states. We can avoid this complication by showing 
that a sequence of transformations can be devised such 
that for each one of them, the initial and final states 
have at least one optimal permutation matrix in com- 
mon. This will allow us to get rid of the independent 
optimization for the initial and final values £^.(i/j). 

Any quantum state with reduced density matrix g 
whose eigenvalues are pi is majorized by a pure state 
with q = (1, . . .), which means that 



\p)=T d _ 1 ---T 2 T 1 \q) 



(B2) 



where the 2Vs are T-transforms [i.e., matrices of the 
form T(i) = (1 — i)l + tW where W is a transposi- 
tion of two particular elements and < t < 1], and 
tk are their respective arguments. We will say that 
T(t) = (1 - t)t + tW is a T-transform of the kind 
W. T-transforms with t > 1/2 can be reduced to T- 
transforms with t < 1/2 by prepending the W permu- 
tation: T{t)W = T(l - t). Thus, Eq. ((B2]) can be cast 
as 



|p) =T d _ 1 W d - 1 ---T 2 W 2 T 1 W 1 \q), 



(B3) 



where W\ , W 2 , . . . are appropriate permutation matrices 
(could be the identity) and now all T-transforms have 
arguments tk < 1/2. 

Moreover, any T-transform with < t < 1/2 can be 
split into an arbitrary number of intermediate transforms 
of the same kind T(t'),T(t"), . . . with < t',t",... < 
1/2. This can be seen by re-paramctcrizing the set of 
T-transforms of any given kind [0 < t < 1/2] by 



T(s) 



1 



-1 



1 



-W, 



(B4) 



where T(0) = 1 and < s < oo. Ob- 
serving that T(si)T(s2) = T(si + s 2 ) shows that 
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with positive sum- 
This allows to de- 



product of T£ k 



nn 



any splitting of s = JA s 
mands can be accomplished, 
compose each Tk into a 
such that the states (Tjr 1 ■ ■ -T^T^! ■ ■ ■T 2 T 1 \q) and 

( T k T k~ 1 ■ ■ ■ T k) T fe-i ' ' ' T 2Ti\q) share at least one opti- 
mal permutation matrix, for all i, k. Thus, finally the 
target state \p) can be written as 



\n)w k 



k=l L \i=l 



\q) 



(B5) 



We are now ready to show that Eq. (|B1I) holds for any 
transformation \p) — > Tj.\p). Let a be a common optimal 
permutation matrix for both states, and A a = <r T Aa, 
and T l k = (1 - t)t + tW, with t < 1/2. The increment in 
£ ^ can be cast as 

= (p|A CT |p)-(p|T(i) T A CT T(Ob) 

= (p| ((1 - (1 - tf )K - t(l - t) [A a W + WK]) \p) 

-t 2 { P \wKw\p) 



>2i(l-t)(p| 



-(VLA CT + A CT VL) 



If)- 



(B6) 



In the last inequality we have exploited the assumption 
that a is the optimal permutation, hence (p|A CT |p) > 
(p\WA a W\p). 
Evaluating expression (|B6[) we find 



(Pi 



A CT — — (VLA CT + A a W) 



b) = 2Ap 2 sin 2 — , (B7) 



where Ap 2 = (p.; —pj) 2 where i, j are the indices swapped 
by W, or the subspace where T\. acts non-trivially, and 
AB = 6 a -i - 6 a -i = ^(o-r 1 - vj 1 )- Indeed, since 



A9 



(B8) 



we have 



AS^{p) > 4£(1 - t)Ap 2 sin 2 



(B9) 



It is then finally straightforward to verify that 4t(l — 
t)Ap 2 = 2^AE L . ' □ 
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